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Using Monte-Carlo methods, we determine the best-fit 
value of the homogeneity parameter a in the Dyer-Roeder 
distance-redshift relation for a variety of redshifts, inhomo- 
geneity models and cosmological parameter values. The re- 
lation between a and the fraction of compact objects, /p, 
is found to be approximately linear. This relation can be 
parametrized with reasonable accuracy for all cases treated 
in this paper by 1 ~ a = a ■ fp, where a ~ 0.6. 

PACS numbers: 98.62. Sb, 98.62.Py, 95.30.Sf 



I. INTRODUCTION 

Assuming the validity of the cosmological principle, 
i.e., that the Universe is isotropic and homogeneous, we 
can apply the Robertson- Walker (RW) metric to the field 
equations of general relativity to derive a relation of cos- 
mological distance measures to redshift for various cos- 
mological parameter values. However, we know that our 
Universe is very far from homogeneous on scales smaller 
than galaxy clusters. It is generally assumed that this 
does not affect the large scale expansion rate of the Uni- 
verse. Still, inhomogeneities will affect measured dis- 
tances through the effect of gravitational lensing. There 
will be different amounts of matter along different lines- 
of-sight, causing different amounts of focusing of the 
light-rays. It is not possible to obtain exact solutions 
of the field equations for general inhomogeneous models, 
thus one is referred to numerical simulations to compute 
gravitational lensing effects on distance measurements. 

Sometimes, we would like to be able to use simpler 
methods to compute at least approximate distances. The 
Dyer-Roeder (DR) distance-redshift relation ||l| assumes 
that the expansion rate of the Universe is governed by 
the total matter density whereas the focusing of light is 
only affected by a fraction a of the total matter density. 
The DR distance thus contains an additional parameter, 
namely the homogeneity-parameter, a. The approxima- 
tion should be fair if a fraction 1 — a of the matter density 
is in very compact objects and the light-ray travels far 
from all matter accumulations, i.e., if one can neglect the 
effect of gravitational lensing. 

In this paper we investigate properties of the DR 
distance-redshift relation by comparing with numerical 



simulations of the distance-redshift relation in inhomo- 
geneous universes, including the effect from gravitational 
lensing. More specifically, we compute the best fit value 
of the homegeneity-parameter a for different cosmolo- 
gies and inhomogeneity models. These values can be 
used, e.g., with some of the publicly available routines 
for computing cosmological distances . 

In an earlier study, Tomita |^ has found the best-fit 
value of a to be one in most cases, with a dispersion in 
a dependent on the cosmological model, the redshift and 
separation of light-rays. 



II. METHOD AND RESULTS 

In Fig.|l], we have plotted the DR angular-diameter dis- 
tance in units of the Hubble-length for different values of 
the homogeneity-parameter a. The larger the value of 
a, the larger the amount of focusing of the light-rays. 
Thus, the angular-diameter distance decreases monoton- 
ically with a, i.e., it increases with the " inhomogeneity- 
parameter", 1 — a. Note also that there is a simple 
relation between the luminosity distance, c^l, and the 
angular-diameter distance cJa, 



dL = (1 + zf ■ dA. 
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FIG. 1. The Dyer-Roeder angular-diameter distance in 
units of the Hubble-length for different values of the homo- 
geneity-parameter a and (Omj^^a) = (0.3,0.7). 

Using the simulation package SNOC Q|, we have used 
light-ray tracing to obtain angular-diameter distances 
with different inhomogeneity models and different values 
of the cosmological parameters. Gravitational lensing ef- 
fects are calculated by integrating the geodesic deviation 
equation through a number of consecutive cells between 
the observer and the source. In each cell, we can spec- 
ify the matter distribution governing the deviation. For 
more details of the method, originally proposed by Holz 
and Wald, see [||J^. In inhomogeneous models, there will 
not be a one-to-one correlation between the redshift and 
the distance since gravitational lensing will cause a dis- 
persion in the Hubble diagram. Figure || shows one sim- 
ulated data set of angular-diameter distances together 
with the DR angular-diameter distance for three differ- 
ent values of a. Using x^-tests, we determine the best-fit 
a-value for each of our simulated data sets. 
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results, see 

We have used three different sets of cosmological pa- 
rameter values; one open with (fiMji^A) — (0.2,0) 
and two flat models with (SlM,f^A) = (0.3,0.7) and 
(r^Mi^^A) = (IjO), respectively. 

In Fig. H to H we present results for (ilM,f^A) — 
(0.3, 0.7) with one homogeneous component and one com- 
ponent in point-masses. Distances are calculated for 
z = 1, z = 2 and a distribution of redshifts, 0.1 < z < 3. 
For all models and redshifts, there is a linear relation 
between the fraction of the total matter density in point- 
masses, /p, and the inhomogeneity-parameter, 1 — a. In 
the simulations using homogeneously distributed matter, 
we have added the condition that 1 — a = for /p = 
when fitting the linear function. 
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FIG. 3. Results for the homogeneity-parameter a for the 
case with one homogeneous component and one component 
in point-masses for 2 = 1 and (JIm, ^a) = (0.3, 0.7). 



FIG. 2. The dispersion in dA for (f2M, Oa) = (0.3, 0.7) with 
a fraction 0.6 of the matter density in galactic dark matter 
halos and a fraction 0.4 in compact objects. 



The case closest to the premises of the derivation of the 
DR distance-redshift relation is the case with one com- 
ponent of the matter density homogeneously distributed 
and another in very compact objects, e.g., point-masses. 

A perhaps more realistic model of our Universe has 
one part of the total matter density in compact objects 
and another part in some realistic galaxy dark matter 
halo model, e.g., the Navarro- Frenk- White (NFW) den- 
sity profile Note that the exact parametrization of 
the galaxy density profile does not significantly affect the 
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FIG. 4. Results for the liomogeneity-parameter ct for the 
case with one homogeneous component and one component 
in point-masses for z = 2 and (SIm, Oa) = (0.3, 0.7). 
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FIG. 6. Results for the homogeneity-parameter a for the 
case with one component with the NFW density profile and 
one in point-masses for 0.1 < z < 3 and (SIm, Oa) = (0.3, 0.7). 
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FIG. 5. Results for the homogeneity-parameter a for the 
case with one homogeneous component and one component 
in point-masses for 0.1 < 2; < 3 and {yim,0,A) = (0.3,0.7). 

In Fig. H to |[ results for the more realistic case with 
one part of the total matter density in point-masses and 
another in dark matter halos parametrized by the NFW 
density profile are presented. Distances are calculated in 
a broad redshift interval 0.1 < z < 3. 
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FIG. 7. Results for the homogeneity-parameter a for the 
case with one component with the NFW density profile and 
one in point-masses for 0.1 < 2 < 3 and (r2M,nA) = (0.2,0). 
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FIG. 8. Results for the homogeneity-parameter ol for the 
case with one component with the NFW density profile and 
one in point-masses for 0.1 < 2 < 3 and (r2M,nA) = (1,0). 




1500 



1000 



1 



n„=0.3,Q,=0,7 



1 .5 



2.5 



3 



FIG. 9. Number of cells traversed in the ray-tracing simu- 
lations for (^M.f^A) = (0.2,0) (top), (^M.f^A) = (1,0) (mid- 
dle) and (^m.^a) = (0.3,0.7) (bottom). 



III. DISCUSSION 

For light-rays passing far from all point-masses, we ex- 
pect the best fit value of a to be given by 1 — a « /p 
where /p is the fraction of the total matter density in 
point-masses. Due to the magnification from lensing we 
expect that on average, 1 — a < /p. It is evident from 
Fig. that the relation between a and /p is indepen- 
dent of redshift. However, from Fig. |^||, we can see that 
there is a weak cosmology dependence. We can under- 
stand this effect from the fact that lensing effects should, 
for a given redshift, be roughly proportional to the dis- 
tance travelled by the light-ray as well as the mass of 
the lenses. In Fig. ^, the number of cells traversed in 
the ray-tracing simulations are presented as a function of 
redshift for the three different cosmologies used in this 
paper. Since the number of cells is proportional to the 
distance travelled, we see that the (r^Mj^^A) = (0.3,0.7) 
cosmology yields the farthest distance, followed by the 
(Om, ^^a) = (0.2, 0) and the (f^M, ^^a) = (1, 0) case. How- 
ever, since the mass in lenses is proportional to J^m, we 
expect that we will have the largest lensing effects for 
{VLyv.^K) = (1,0) followed by (Slyv.'^K) = (0.3,0.7) and 
(r^M, ^^a) = (0.2, 0), as is confirmed by Fig. 



We can estimate the dispersion in a from Fig. ^ to 
be fJa ~ 0.2. However, since we are measuring dis- 
tances, the dispersion in distance at a given redshift is 
perhaps of greater interest. In Fig. the dispersion in 
dA for /p = 0.2 and (SlyiMh) = (0.3,0.7) for z = 0.5, 
z — I and z = 1.5 is shown. The zero- value corresponds 
to the value one would obtain in a homogeneous uni- 
verse, the so called Robertson- Walker angular-diameter 
distance, dA.RW- The dispersion in angular-diameter dis- 
tance (rms) corresponds to 5.9/h Mpc for z = 0.5, 22/h 
Mpc for z = 1 and 43//i Mpc for z = 1.5. It is also evi- 
dent from Fig. |l^ that the RW value is within the rms- 
dispersion in all three cases. However, since the mean 
value is displaced from zero, the use of the RW-distance 
will cause systematic errors of the order 3.5/ h Mpc for 
z = 0.5, 13/h Mpc for z = 1 and 24/ft. Mpc for z = 1.5. 
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FIG. 10. The dispersion in dA for the NFW case with 
/p = and {^mMa) = (0.3,0.7) for z = 0.5, « = 1 and 
z = 1.5. 



IV. SUMMARY 



We have investigated the best-fit values of the 
homogeneity-parameter a in the Dyer-Roeder distance- 
redshift relation. For a variety of inhomogeneity models, 
redshifts and cosmological parameter values, the relation 
between a and the fraction of compact objects, /p, is ap- 
proximately linear and we can parametrize the relation 
by 

1 - a = a • /p, 

where a ~ 0.6. We expect the dispersion in the angular- 
diameter distance to be 5.9/h Mpc for z — 0.5, 22/h 
Mpc for z = 1 and 43/h Mpc for z = 1.5 in doing this 
approximation. 
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